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Abstract 



^ . In this paper we present a method to improve the description of 0+1 dimensional boost 

invariant dissipative dynamics in the presence of large momentum-space anisotropics. We 
do this by reorganizing the canonical hydrodynamic expansion of the distribution func- 
tion around a momentum-space anisotropic ansatz rather than an isotropic equilibrium one. 
^ ■ At leading order the result obtained is two coupled ordinary differential equations for the 
C^. momentum-space anisotropy and typical momentum of the degrees of freedom. We show 
that this framework can reproduce both the ideal hydrodynamic and free streaming lim- 
^ \ its. Additionally, we demonstrate that when linearized the differential equations reduce to 
Q^ \ 2nd order Israel-Stewart viscous hydrodynamics. Finally, we make quantitative conipar- 
^ ■ isons of the evolution of the pressure anisotropy within our approach and 2nd order viscous 
Q . hydrodynamics in both the strong and weak coupling limits. 

P^ \ Keywords: Boost Invariant Dynamics, Anisotropic Plasma, Non-equilibrium Evolution, 
Q \ Viscous Hydrodynamics. 

---J . 1. Introduction 

X 

Jj \ The remarkable success of relativistic hydrodynamics to describe the anisotropic flow of 

matter created during non-central heavy ion collisions [ij, |2|, y, Ij] has sparked a great deal 
interest in the systematic derivation and application of viscous hydrodynamical equations 
for relativistic systems p, Ig, |7|, |8|, lol, |lO|. Here we demonstrate a method for deriving 
hydrodynamic-like equations for systems which have large momentum-space anisotropics. 
Such large momentum-space anisotropics can arise naturally due to the rapid longitudinal 
expansion of the matter created during relativistic heavy ion collisions. These anisotropics 
can become so large that the shear becomes as large as the isotropic pressure, signaling the 
breakdown of the expansion of the energy-momentum tensor in terms of shear corrections. 
In fact, in 2nd order viscous hydrodynamics large momentum-space anisotropics can even 



cause the longitudinal pressure predicted to become negative [UJ 
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Here we present a way to circumvent such problems and derive evolution equations which 
can describe the dynamics of systems with potentially large momentum-space anisotropies. 
The method is to change the usual hydrodynamic expansion of the one-particle distribution 
functions 

/(r,x,p) = /eq(|p|,T(r)) + ,5/i + 5/2 + --- , (1) 

which is expansion about an isotropic equilibrium state /eq(|p|, T(r)), to one in which the 
expansion point itself can contain momentum-space anisotropies 

/(^, X, p) = /aniso(P,Phard(r), ^(r)) + 6 f[ + 6 f2 ^ , (2) 

where in both equations r is proper time. In Eq. ([2]) ^ is a parameter which measures 
the amount of momentum-space anisotropy and phard is a non-equilibrium momentum scale 
which can be identified with the temperature of the system only in the limit of isotropic 
equilibrium. We will introduce an ansatz for /aniso which approximates the system via 
ellipsoidal equal occupation number surfaces as opposed to the spherical ones associated 



with the isotropic equilibrium distribution function [12j . For highly anisotropic systems one 
immediately expects that the corrections 5/^ will have smaller magnitude than the isotropic 
corrections 6fn due to the fact that momentum-space anisotropy is built into the leading 
order of the expansion in the reorganized approach. One can use different prescriptions for 
an anisotropic basis. Here we have in mind the use of spheroidal harmonics \l^ for which 
our ellipsoidal ansatz is the leading order term. 

Using this as a starting point, we will derive leading order 0+1 dimensional boost- 
invariant evolution equations for /aniso by taking moments of the Boltzmann equation]^ The 
result will be two coupled ordinary differential equations which describe the evolution of the 
system's momentum-space anisotropy and typical hard momentum scale. We demonstrate 
that at leading order the new expansion can equally well describe both the ideal hydrody- 
namic and free streaming limits within a unified framework. In addition, we demonstrate 
analytically that when the associated coupled nonlinear differential equations are expanded 
to leading order in the anisotropy parameter, they reduce identically to the 2nd order viscous 



hydrodynamic equations of Israel and Stewart [Ij, |15|, |16|, ll7[. In the general case we solve 
the coupled nonlinear equations numerically allowing us to describe both the early-time 
high-anisotropy dynamics and late-time near-equilibrium dynamics using a single formalism 
which gives us as output the time-evolving distribution function. This distribution function 
can then be used as input for other calculations. 

2. Deriving boost invariant dynamics from transport theory 

We consider a highly populated system formed by massless particles. We describe 
this system by means of the Boltzmann equation for the one-particle distribution func- 
tion /(t,p,x) in the lab frame. Additionally, we assume that the system is boost invariant 
and expands only along the longitudinal (beam line) axis. Boost invariance implies that 



"'^By leading order, we mean neglecting the corrections from the Sf^ terms in Eq. ([2|). 
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that the longitudinal velocity of the system is constant along lines of constant spatial rapid- 
ity such that Vz = z/t. The assumption that the system only expands in the longitudinal 
direction is reliable if the transverse size of the system is sufficiently large (~ 1.2 A^/^ for 
central nucleus-nucleus collisions) that the effects of transverse dynamics can be ignored to 
first approximation. If this is the case then one can assume a homogeneous distribution in 
the transverse directions and set v^ = Vy = Q. With these assumptions one can show that 
the dynamics reduces to 0+1 dimensional evolution, i.e., only derivatives with respect to 
proper time remain. Here we will show how this occurs by starting in lab coordinates and 
transforming the Boltzmann equation to comoving coordinates. Under the assumption of 
homogeneity in the transverse direction one has in lab coordinates 

p'dtf{t, z, p) + p^dj{t, z, p) = -C[/(t, z, p)] , (3) 

where C[/(t, z, p)] is the collisional kernel which is a functional of the distribution function. 
Its functional form depends on the type of interactions between the particles and, in general, 
is difficult to evaluate exactly even for simple theories. To make analytic progress in this 
paper we will study the case that the collisional kernel is given by the relaxation time 
approximation 

C[/(t, z, p)] = pX r [fit, z, p) - /eq(t, z, IpI, T(r))] , (4) 

where /cq is the local equilibrium distribution function and F is the relaxation rate. The 
temperature in the local rest frame T{t) is determined dynamically by requiring ^eql'T") = 
^non-eqij), that guarantees energy conservation within this approximation [l8| and T(r) is 
assumed to only depend on the proper time r = \Jt'^ — z^ in accordance with the assumption 
of boost invariance. Although here we consider the relaxation time approximation, we note 
that the general method can be applied to collisional kernels derived directly from quantum 
field theoretical methods. We postpone such studies to future publications. 

2.1. Ansatz for the distribution function and equation of state 

In order to have a tractable approach that will allow us to describe systems that may be- 
come highly anisotropic in momentum space we consider a general ansatz for our anisotropic 
distribution function which was first introduced by Romatschke and Strickland (RS) |l2| . 
Within this approach the leading order anisotropic distribution function in the local rest 
frame can be obtained from an arbitrary isotropic distribution function (/iso) by squeezing 
[C, > 0) or stretching (^ < 0) /iso along one direction in momentum space 

/(r,x,p) = /«5(p,e(r),Phard(r)) = /iso([p' + e(r)(p- n)VpL.dW) , (5) 

where Phard is related to the average momentum in the partonic distribution function, fi is 
the direction of the anisotropyp and — l<^<ooisa parameter that refiects the strength 
and type of anisotropy. The anisotropy parameter ^ is related to the average longitudinal 



and transverse momentum of the constituents via the relation [19|, |20 



'^m-'- <«' 



^Hereafter, we use n = e^, where e^ is a unit vector along the longitudinal direction. 
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The system is locally isotropic when ^ = but this does not imply that it is in local thermal 
equilibrium unless /iso is an equilibrium distribution function. 

The energy-momentum tensor in the local rest (comoving) frame with coordinates x^ = 
(r = Vt"^ — z^,x,y,c; = arctanh(2;/t)) is given by T^" = (27r)~^ J(i^p/p°p^p^/(r, x, p). In 
the comoving frame the energy-momentum tensor is diagonal and using the RS ansatz i^ 
its components are 



21 



^(Phard, 0=T^^ = 2 ( YTJ ^ 7f ) ^i«°(Phard) , (7a) 

= 7^(0 ^iso(Phard) , 

^T(Pha.d, = 1 (T^^ + Tyy) = ^ (^i±i^^lMi)^ vuv^^r,) , (rb) 

= 7^T(0^iso(Phard), 
^L(Phard,0 = -^c = 7 TT—. 7^iso(Phard) , (7c) 



= '^L(0^iso(Phard), 

where PisolPhard) and £^iso(Phard) are the isotropic pressure and energy density, respectively. 
In general, one cannot identify immediately Piso(Phard) and ^iso(Phard) with their equilibrium 
counterparts, unless one implements the Landau matching conditions |2l|.|f| For a conformal 
system the tracelessness of the stress-energy tensor T/^=0 implies £ = 2Vt + Vl- This 
condition is satisfied by Eqs. ([7]) for any anisotropic distribution function ([5]) since for an 
isotropic conformal state £^iso = SPiso- 

In addition to the components of the energy-momentum tensor, one can also calculate 
the entropy density in the comoving frame following the standard definition from kinetic 



theory [23| . For the case of the RS ansatz (^ one obtains |21 

S = -/^/(r,x,p){log[/(r,x,p)]-l} 



SisoiPl 



hard J 



where iSiso(Phard) is the isotropic entropy density. 

An additional property of the RS ansatz (|5]) is that in the case of 0+1 dimensional boost 
invariant expansion, it is possible to relate the anisotropy parameter ^ and the shear used in 
viscous hydrodynamical treatments by matching the values of the longitudinal pressure Pl- 
To see this, consider the longitudinal pressure for 0+1 dimensional viscous hydrodynamics 
with an ideal equation of state 



■^In general the Landau matching conditions are a way to connect non-equilibrium evolution to near- 
equilibrium hydrodynamic evolution. In practice, this cor resp onds to matching the components of the 
energy-momentum tensor and matching collective velocities [2, 
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where 11 = 11^ is the ^<j component of the fluid shear tensor. The local temperature T{t) and 
hard momentum scale phard contained in the RS ansatz ([5j) are related through the Landau 
matching condition for the energy density, Scq^r) = Sn^-eq^r), which implies that at any 
time r we can relate T and Phard via T = TZ^^^iOPhard (2l|. By using this relation, Eq. (17c|) 
can be rewritten as 

Vdr) = 7^L(e(r))P^(Phard(r)) = ^#;^P;^(T(r)) 



geq(T(r)) 7^L(e(r)) 
3 7^(e(r)) 

Equating Eqs. ([9]) and (ITO|) we find 



(10) 



u(r) ^ ^a^l)) 



1 



^L(e(r)) 



:iii 



which is valid to all orders in C,- In the case of small departures from equilibrium the system 



is very nearly isotropic with 11 ^ £^ and one can use a small ^ expansion to obtain [21 
where £eq is the equilibrium energy density. If one uses the Navier- Stokes value of the shear 



tensor IIns = 4?7/(3r) and the ideal equation of state in the last relation, we obtain [24 



^Ns = ^1 + 0(n^s) , (13) 

where r] is the shear viscosity and S is the equilibrium entropy of the system. Therefore, 
deviations from equilibrium which are canonically inferred from the evolution of the shear 
viscous tensor can be also studied using the RS distribution function (|5]) by analyzing the 
evolution of the anisotropy parameter ^. 

We point out that Phard and ^ are time-dependent variables whose evolution equations 
are unknown a priori. In this work, we derive and solve such equations by taking the first 
two moments of the Boltzmann equation within the relaxation time approximation. This 
procedure allows us to study their time evolution and obtain information about the tran- 
sition between early-time non-equilibrium dynamics and late-time viscous hydrodynamical 
behavior. 

2.2. Boltzmann equation in comoving coordinations 
It is convenient to switch to comoving coordinates 

t = r cosh <;- , 

z = rsinh<;", (14) 



and introduce momentum-space rapidity 

p° = Pt cosh y, 
P^ = Pt sinh y . 



In this coordinate system the metric Qfj^^, = diag 
is /Rs(p,^,Phard) = /iso (pt[1 + (1 + sinh^(y 
Boltzmann equation ([3]) can be written as 



(15) 

—1, — T^) and the RS ansatz (^ 
^)]/Phard)- After changing variables, the 



.dpT sinh(y-^) d\ 

PTCOSh{y-q)—^ — 1 /Rs(P,4,Phard) 

-TpT cosh{y - ^)[/Rs(p,^,Phard) - /cq(p, T(r))] 



(16) 



Note that in the above equation one should evaluate the last term on the right side using 
the constraint T{t) = R^^'^{^)pha.rd which is required by energy conservation la, |21 . 



The usual way to study non-equilibrium deviations within a kinetic theory approach is 



to take moments of the Boltzmann equation J25|. For instance, by using the 14 moments 



ansatz [IJ, ll5|] it is possible to determine the equation of motion for the shear viscous 
tensor by taking the 2nd moment of the Boltzmann equation. In this work, we just need 
to determine the equations of motion for ^ and phard by considering the zeroth and first 
moments of the Boltzmann equation flT6|l . We demonstrate this in the next two sections. 



2.2.1. Zeroth moment of the Boltzmann equation 

The zeroth moment of the Boltzmann equation with respect to the single particle energy 
reduces to dfj^N^ on the left hand side 23|; however, in the relaxation time approximation the 
zeroth moment of the collisional kernel on the right hand side is non- vanishing. Since in any 
realistic non-equilibrium quantum field theory there are generically number non- conserving 
processes we consider this to be a useful feature of the relaxation time approximation, p 

To derive the first equation, let us first rewrite the left hand side (LHS) of Eq. ( 1T6|) 
explicitly using the RS ansatz 



LES = Ed^Us{co) 



' pI 

Phard 



dr^- 



2(1+0 



2u 

Phard 



drPh, 



ard 



(17) 



where u = [p^ + (1 + ^)pl] /Phard ^^^ ^^ have rewritten the result in terms of the local 
energy and longitudinal momentum defined in the fluid cell via E = pt cosh.{y — i;) and 
Pz = pTsinh.{y — q), respectively. By using this expression together with the right hand side 



'^We note that if one were considering a number-conserving theory, one can enforce 9^iV'' = in the 
relaxation time approximation by introducing a particle fugacity into the ansatz for the one particle distri- 
bution function. One would then need to take one further moment of the Boltzmann equation to determine 
the dynamical equations. 



of Eq. (TT6l) we can take the zeroth moment by integrating both sides with the momentum- 
space measure {2tt)^^ J d^p/{2p^). This results in 



(pld^fRsiuj)) 



Phard 



dr^ 



2(1 +0\ 2{ud^fKsico)) 



-r 



drPhard 
Phard 

/Rs(P,'^,Phard)) - (/cq(p, T(r))) 



where (a) = (27r) ^ Jd^pa{p). These integrals can be calculated analytically using the 
general form of the RS ansatz i^ 



(pld^Msiuj)) 
(ujduifRsiuj)) 

(/Rs(P,C,Phard)) 



Phard 



1 

"2 + 0=^/2 

"J '^isolPhardj 

"2 (1 + 0^/2 

"iso l^Phard / 



riisoiP: 



hard/ 



(1 + ^)1/2' 
(/eq(p,T(r) = 7^1/'(0Phard)> = 7^=^/'(0 ^iso(Phard) , 

where in the last line we have used the Landau matching condition T{t) 
After replacing these expressions in Eq. fllSp . we obtain 
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^ + ^ r Phard 



2r 



n'/\oVi + ^ 



(19a) 

(19b) 

(19c) 
(19d) 

7^^/'(0Phard. 

(20) 



which is our first coupled differential equation for ^ and phard- 



2.2.2. First moment of the Boltzmann equation 

In this section we derive the second coupled differential equation for ^ and phard- Taking 
the first moment of the Boltzmann equation is equivalent to energy-momentum conservation 
dfj,Tf^'^ = |23|. We enforce energy conservation by requiring that £eq{T) = ^non-eq('7") which 
results in T{t) = R^^^{^)phard- Momentum conservation is guaranteed for any distribution 
function which is parity-symmetric in momentum space. In the comoving frame, energy 
conservation implies that the equation of motion for the energy density is 

dSjr) _ S{t)+V,{t) ^ 
dr T 

where S and Vl are the energy density and longitudinal pressure. Using the components of 
the energy momentum tensor given by the RS ansatz (Eqs. (I7ap and (17c]) ) we obtain 



Phard 



where 



mi 

7^(0 

7^'(0 = d^nio 



ard 



1 






— 1 



arctan-\/^\ 



Eq. (122|) is the second coupled ordinary differential equation for ^ and phard- 
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(22) 



(23) 



2. 3. Combined Differential Equations 

It is possible to further simplify Eqs. (120]) and ( 12^ to obtain 






(24a) 



1+4 Phard 27^(4) + 3(1 + 4)^ (0 

In the results section we show that the solution of these equations allows us to describe 
the transition between early-time longitudinal free streaming and late time hydrodynamical 
behavior. Note that one should expect that the solutions of these equations reduce to 
either free streaming expansion or ideal hydrodynamical behavior when F — )■ or F — )■ oo, 
respectively. In the next section we show this by taking the asymptotic limits of these 
equations. 

2.4- Asymptotic limits 

One way of checking self-consistency in our approach is to see if it is possible to obtain 
from Eqs. (I24p the two extreme cases of expansion: (1) a system which undergoes 0+1 di- 
mensional longitudinal free streaming expansion and (2) a system which is "instantaneously" 
thermal and isotropic and undergoes 0+1 dimensional ideal hydrodynamical expansion. In 
this section we show that Eqs. fl2^ contain both of these limits naturally. 



2.4- 1- Free Streaming Limit 

When the particles do not interact the interaction rate F vanishes exactly. In order to 
show that the free streaming limit is recovered in our approach, we set F = in Eqs. (1241) . 
One obtains immediately 

dr^ = -{l + 0, (25a) 

r 

drPh.rd = . (25b) 

The solutions of these equations are 

^{r) = il + ^o)(jr] -1, (26a) 

Phard = PO , (26b) 

where ^^o and po are the initial values of ^ and phard at r = tq. This corresponds precisely to 



the free streaming solution to the Boltzmann equation (18|, |20|, l2l|, |26 . 



2.4-2. Ideal hydrodynamical expansion 

In the case that the system is isotropic at all times we have ^(r) = and therefore 
dr^ = 0. This occurs in the limit that F — ;■ oo. From Eq. ( 120|) . we can expand the right 
hand side 



lim2F 



R'/\ov^+^\ = 2F[i-(i+o(a)] , 

= 0, (27) 

so that in this limit Eq. fl2CT]) gives us 

^rPhard = "— • (28) 



"-' 7 /^ 11(11 U (-, 

Phard 3r 

We can also consider the same limit in Eq. (E2|) . First, we expand the right hand side 
so that Eq. (l22l) reduces to 

^rPhard = "^ • (30) 

Eqs. fl28|) and fl30l) are exactly the same and the solution is Phard(T) = iOhard(rn) (rn/r)^/'^ 



which corresponds to the solution for the ideal 0+1 boost invariant expansion [27 



2. 5. Positivity of entropy divergence 

One can prove analytically that the solutions to Eqs. (I24p have positive entropy diver- 
gence. Formally the requirement is stated as d^S'^ > where S^ is the entropy current. 
In the 0+1 dimensional case this simplifies to the requirement drirS) > where S is the 
anistropic entropy defined in Eq. (|8]). Using Eq. ([8]) we first use the chain rule to obtain 

S 2(1 + g Phard 

Using Eqs. f l2^ and simplifying we find 



tT 



n'/\OV^+l-^ ■ (32) 



Since 7?.^/^(^)vT+~^ > 1 for all ^ in — 1 < ^ < oo and all other quantities are positive, we 
therefore always have drirS) > 0. Note that in the ideal hydrodynamical limit (^ = 0) and 
the free streaming limit (F = 0) the lower bound of the inequality is reached with drirS) = 
as is expected. 



^We point out that in the limit ^ H> 0, the local equilibrium temperature T{t) is exactly the same as the 
typical momentum of the particles in the system phard since 7?.(0) = 1. 
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2. 6. 2nd order viscous hydrodynamics from RS ansatz 

In Sect. 12. II we showed that in the near-equihbrium hmit one can relate the viscous shear 
tensor 11 with the anisotropy parameter ^ (Eq. flT2l) ). As a consequence, in our approach, 
the evolution of dissipative corrections to the ideal fluid is equivalent to the evolution of the 
anisotropy parameter ^. To see this more in detail, in this section we derive the equations 
of motion for 2nd order 0+1 boost invariant viscous hydrodynamics from the RS ansatz i^ 
by linearizing the previously obtained coupled nonlinear differential equations. 

We begin by using Eq. fl2Tl) for the energy density, Eq. (TTOl) for the matching between ^ 
and n, and assume an ideal equation of state|j The result is 



which upon using (ITTI) gives 



^'^ = -^i^^3#yj- i^^*) 



a^f = _£±Zh.5. (34) 

T T 



This expression is the well-known equation for the energy density when the system has shear 



viscous corrections in a 0+1 dimensional boost-invariant expansion [IJ, [151, |16 . 

To obtain the equation of motion for the shear viscous tensor we use Eq. (I20p derived 
in Sect. 12.2.11 First, let us rewrite some terms which appear in this equation by using the 
linear relation between ^ and 11 

where S is understood to be the isotropic equilibrium energy density in this section. Next, 
by implementing the Landau matching condition for the local temperature and phard we 
obtain 

drP..r,_lfd.£ ^'^%eV (36) 



Phard 4 V ^ 7^(0 

Plugging the last two expressions into Eq. f l20|) . using Eqs. f lT2|) and flMl) . and expanding to 
the lowest non-vanishing order in ^ one obtains 

a^n + i5_15£ = _irn. (37) 

3 r 45 r 2 ^ ^ 

To finally connect with second order viscous hydro, we identifju 

r = — , (38a) 

r. = ^^, (38b) 



^During the derivation presented in this section, it is necessary to implement the Landau matching 
condition which relates the local temperature and the average momentum of the particles phard, i-e- T{t) — 
7^l/4(e(r)Kard(r). 

^This is analogous to the identifications of the unknown coefficients /Sq, /3i and /32 which appear in the 
14 Grad's ansatz [H M, M, [13, S Hi] ■ 
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Figure 1: Hard momentum scale phard as a function of proper time for five different values of 77: the ideal 
hydrodynamic limit (fj = 0), 77 = 0.1, 77 = 1, ^ = 10, and the free streaming limit {fj = cx)). 



SO Eq. fl37|) can be written finally as 

drll = 



n 



4 T] 



4n 

37' 



(39) 



where t^^ and 1] are the shear relaxation time and shear viscosity, respectively. This last 
expression is p recisely the 0+1 dimensional 2nd order viscous hydrodynamical equation for 
n 2A, ll5|, ll6|. Therefore, the physics of 2nd order viscous hydrodynamics can be recovered 
from our ansatz in the small ^ limit. 



3. Results 

In this section we present the results of numerically integrating the coupled nonlinear 
differential equations using the relation between the relaxation rate F and the ratio of shear 
viscosity to equilibrium entropy fj = rj/S which results from Eqs. ( 138|1 . namely 



2T(r) 27^l/4(OPhard 



5?7 



5r] 



(40) 



which should be used for T in Eqs. fl24j) . In all plots shown in this paper we will use the 
following initial conditions: ro= 0.2 fm/c, Phard(T = tq) = 350 MeV, and ^(r = tq) = 0. 
Note, however, that one can choose whatever initial conditions one likes in order to integrate 
Eqs. dMD. 

In Fig. [T] we show the solutions obtained for the hard momentum scale Phard as a function 
of proper time for five different values of f/: the ideal hydrodynamic limit {fj = 0), f] = 0.1, 
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Figure 2: Anisotropy parameter ^ as a function of proper time for three different values of 77: 77 — 0.1, 
77 = 1, and 77 — 10. 



fj = 1, fj = 10, and the free streaming limit (f/ = 00). As can be seen from this figure, 
the solutions to Eqs. fl2^ make a smooth transition from ideal hydrodynamic behavior to 
longitudinal free streaming as f/ goes from to cxd. For finite values of f] one sees that the 
solutions smoothly transition between early-time longitudinal free streaming and late-time 
viscous hydrodynamical behavior. 

In Fig.|2]we show the corresponding solutions for the anisotropy parameter ^ as a function 
of proper time for three different values of r;: fj = 0.1, fj = 1, and f/ = 10. In the ideal 
hydrodynamic limit one has ^ = and in the free streaming limit one has ^ = (l+Co)(''"/'^o)^~ 
1. The solutions presented in Fig. [2] show that at very early times the behavior of the system 
only depends weakly on fj due to the rapid longitudinal expansion of the system. In addition, 
in all cases, the anisotropy has a maximum which is obtained at later proper times as fj is 
increased. 

Knowledge of the proper-time dependence of Phard and C, allows one to calculate all of 
the various thermodynamic functions defined in Eqs. (17al) . (l7b|) . (fTcj) . and ([8]). For sake 
of brevity we do not show each of these individually but we note that the energy density 
smoothly transitions from early-time free streaming to late-time hydrodynamical expansion 
as one would expect. A quantity that is sensitive to momentum-space anisotropics is the 
ratio of the longitudinal to transverse pressures, Vl/Vt- In the ideal hydrodynamic limit 
this is identically one and in the free streaming limit one obtains Vl/Vt oc 1/t^ at late times. 
In Fig.[3]we show Vl/Vt for three different values off/: fj = 0.1, fj = 1, and fj = 10. One sees 
from this figure that as fj is increased, the longitudinal pressure is decreased; however, at 
late times the system restores isotropy in momentum space for any finite fj. It is important 
to note that even for the case fj = 10 the longitudinal pressure of the system remains positive 
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Figure 3: Ratio of longitudinal [Eq. ([7c)) ] and transverse [Eq. (|7b[) ] pressures as a function of proper time 
for three different values of ^: f\ — 0.1, 77 = 1, and f\ — 10. 



definite. 

Having obtained the general numerical solutions to Eqs. fl24p we can now make a quan- 
titative comparison between solution to the fully nonlinear differential equations (12^ and 
the linearized 2nd order viscous hydrodynamical equations given by Eqs. f lM|) and fl39|) . In 
Figs, m and |5] we compare the two approaches for typical strong and weak coupling values 
oi f] = l/(47r) and f] = 10/(47r), respectively. As one can see from Fig. IHfor f/ = l/(47r) 
the difference between the RS ansatz solutions and 2nd order viscous hydro is at maximum 
approximately 15%. However, in the weak coupling case of f/ = 10/(47r) one finds that 
the correction can be greater than 100%. More importantly, one finds that in the weak 
coupling case, for the chosen initial conditions, the ratio Vl/Vt obtained from 2nd order 
viscous hydrodynamics becomes negative for 0.65 fm/c <t< 5.9 fm/c whereas the solution 
obtained using the RS ansatz is positive definite during the entire evolution. Quantitatively 
one should note that even in the strong coupling case large momentum-space anisotropics 
can be developed with TumiVL/VT) ~ 0.45 even for the exact solution. The precise num- 
bers, of course, depend on the choice of the initial conditions and proper time at which one 
begins integrating the differential equations. When cast into dimensionless form the relevant 
quantity is TqT with larger values of the dimensionless number corresponding to generation 
of smaller momentum-space anisotropics. 



4. Conclusions and Outlook 

In this paper we have shown that by starting from an ansatz which incorporates momentum- 
space anisotropics from the beginning one can obtain solutions which smoothly connect the 
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Figure 4: Ratio of longitudinal [Eq. ([7c)) ] and transverse [Eq. (|7b[) ] pressures as a function of proper time 
for fj = l/(47r). Solid line is numerical solution to Eqs. (1^^ and dashed line is numerical solution of 2nd 
order viscous hydrodynamics given by equations (|34p and (j39l) . 



ideal hydro dynamic and free streaming limits by varying fj from zero to infinity. We demon- 
strated that these limits are achieved analytically and numerically. In addition, we showed 
that when the equations are linearized the resulting linear differential equations reduce to 
the Israel-Stewart equations for viscous hydrodynamical evolution. We then made a detailed 
comparison of the prediction for the ratio of longitudinal to transverse pressure predicted by 
the RS ansatz and the Israel- Stewart formalism. We found that for typical strong coupling 
values of f] the correction from 2nd order viscous hydrodynamical evolution was on the order 
of 15% and for a typical weak coupling value of f] that the correction was sizable and could 
be greater than 100%. 

The conclusions reached here are based on one particular set of initial conditions, how- 
ever, the pattern observed here is generic. At LHC energies the initial time of hydro evo- 
lution is expected to be less based on perturbative estimates of the thermalization time; 
however, the initial temperature is also expected to be higher. When cast into dimension- 
less form what matters is the dimensionless combination tqTq for determining the magnitude 
of expected longitudinal momentum-space anisotropies. For the RHIC-like initial conditions 
presented here this combination is tqTq = 0.35 whereas for LHC-like initial conditions one 
obtains TqTq = 0.43 assuming tq = 0.1 fm/c and Tg = 850 MeV. Therefore, one expects 
slightly smaller momentum-space anisotropies to be generated for LHC initial conditions. 

The method shown here can be used as input to calculate the dependence of high energy 
observables on the time evolution of the anisotropy parameter and hard momentum scale. 



There are now calculations of the anisotropic photon rate [301], dilepton rate [19|, |20|, |31 



and quarkonium binding energies J32|, |33|, l3J] using the RS ansatz. It would be interesting 
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Figure 5: Ratio of longitudinal [Eq. ([7c)) ] and transverse [Eq. (|7b[) ] pressures as a function of proper time 
for fj = 10/(47r). Solid line is numerical solution to Eqs. ([M)) and dashed line is numerical solution of 2nd 
order viscous hydrodynamics given by equations ([M)) and (|39p . 

to apply the differential equations derived here to these phenomenological applications. 

In addition, looking forward one can also extend the ansatz used here to allow for the 
inclusion of particle fugacities. This can be done by evaluating another moment of the 
Boltzmann equation; however, one will need to take care with the time evolution of the 
number density in this case if one wants to obtain solutions which reach chemical equilibrium 
at late times |3_5l]. We postpone this development to a future paper. 



Note Added 

During the final preparation of this manuscript another paper presenting a similar idea 
appeared by Florkowski and Ryblewski 36| . However, they obtain two differential equations 
by requiring energy conservation and introducing an ansatz for an entropy source whereas 
we follow the method of taking moments of the Boltzmann equation directly. 
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